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Abstract 



We show how the MHV diagram description of Yang- Mills theories can be used to study 
non-supersymmetric loop amplitudes. In particular, we derive a compact expression for 
the cut-constructible part of the general one-loop MHV multi-gluon scattering amplitude 
in pure Yang-Mills theory. We show that in special cases this expression reduces to 
known amplitudes - the amplitude with adjacent negative- helicity gluons, and the five 
gluon non-adjacent amplitude. Finally, we briefiy discuss the twistor space interpretation 
of our result. 
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1 Introduction 



Since Witten's discovery that topological string theory on super twistor space provides a 
description of A'' =4 super Yang- Mills (SYM) [1], considerable progress has been made 
using twistor-inspired methods to study Yang- Mills theories. An important factor in this 
has been the proposal to use maximally helicity violating (MHV) diagrams, built using 
MHV amplitudes as vertices, in order to derive amplitudes [2]. The MHV diagram con- 
struction has the appealing feature that twistor space localisation is built in. The method 
also carries the great practical advantage of tremendously simplifying the calculation of 
amplitudes. It was quickly confirmed that MHV diagrams gave correct results at tree 
level (see [3] for a review of the field up to August 2004). The prognosis at one- loop was 
initially poor, however, as general arguments indicated that it would be impossible to ig- 
nore conformal supcrgravity fields propagating in the loops [4]. Separately to this, initial 
explorations of the differential equations satisfied by known one-loop amplitudes appeared 
to suggest unexpected complications in their twistor-space localisation properties [5]. 

A direct derivation from MHV diagrams of the one-loop MHV scattering amplitudes in 
A/" = 4 SYM was presented in [6]. The question of the localisation of amplitudes was then 
revisited, and the complications previously found were seen to be due to the appearance 
of additional inhomogeneous terms [7] in the differential equations obeyed by one-loop 
amplitudes. Taking into account the corrections coming from this, one finds that indeed 
the one-loop MHV amplitudes in A/" = 4 SYM localise on pairs of lines in twistor space [8], 
as the direct construction of [6] suggests. This encourages one to conjecture that the whole 
quantum theory of A" = 4 SYM possesses simple twistor space localisation properties. 
By studying the differential equations satisfied by the unitarity cuts of amplitudes, the 
coefficients of the box functions in the next to MHV (NMHV) amplitudes in A/'=4 SYM 
have recently been shown to localise on planes in twistor space [9] . A direct MHV diagram 
construction of these amplitudes has not yet been given however. 

The study of the analytic properties of amplitudes, using the twistor-inspired ap- 
proach, has since been found useful in the general analysis of one-loop amplitudes in 
A/' = 4 Yang- Mills, with a recent derivation of the ( h + + +) one-loop NMHV am- 
plitude [10, 11]. This coincides with one case of the general one-loop seven-gluon NMHV 
amplitude which was also found recently [12] using the cut-constructibility approach. For 
M— 1 theories, the twistor space structure of one-loop amplitudes was studied in [5, 13, 14] 
and it was found that the holomorphic anomaly of unitarity cuts [7] leads to differential 
equations [13], in contrast to algebraic equations for A/" = 4 [10], obeyed by the one-loop 
amplitudes. 

MHV diagrams provide a well-defined prescription for the direct derivation of ampli- 
tudes. It is natural to ask whether the MHV diagram construction of the one-loop N'—A 
MHV amplitudes of [6] can be generalised in other directions - in particular to theories 
with less supersymmetry. This has been confirmed in recent work [15, 16], where the 
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MHV diagram method was shown to correctly reproduce the known MHV amphtudes 
for the M = 1 chiral multiplet. This result implies that one-loop MHV amplitudes for 
all supersymmetric gauge theories can be derived from MHV diagrams, and hence have 
simple localisation properties in twistor space. 

The close relationship between the MHV diagram construction and unitarity-based 
methods [17], first seen in [6], and the success in applying this method to the N =1 case, 
encourages the belief that all cut-constructible amplitudes may be amenable to this new 
approach. It is also of great importance to explore whether MHV diagrams can be used 
at loop level in non-supersymmetric theories.^ These motivations lead one to consider 
the one-loop MHV amplitudes in pure Yang-Mills theory. These amplitudes consist of 
terms containing cuts, which we call the cut-constructible part of the amplitude, plus 
additional rational terms. The amplitudes are of great interest, since they are an example 
of one- loop n-point scattering amplitudes in QCD, where all external particles and the 
particle running in the loop are gluons, and they can be decomposed as 

^glue j^=i ^ ^Ar=l, chiral _|_ ^scalar ^-|^ 

The first term describes the contribution of an A/" = 4 SYM multiplet to the amplitude. 
The second is —4 times the contribution of an A/'= 1 chiral multiplet, and the third is a 
non-supersymmetric amplitude with only complex scalars propagating in the loop. In this 
paper we focus on the calculation of the final contribution since the other two are known. 
A similar supersymmetric decomposition exists for one-loop gluon scattering amplitudes 
with massless quarks or adjoint fermions running in the loop. 

The one-loop MHV amplitude in pure Yang-Mills is known only for two special cases 
- when the two negative- helicity external gluons are adjacent, the cut-constructible part 
is known [19]; and, in the five-gluon case, the full amplitude, including rational parts, 
has been calculated for arbitrary helicity configurations in [20] . In this paper, we will use 
MHV diagrams to derive a compact expression for the cut-constructible part of the general 
one-loop MHV multi-gluon amplitude when there are scalar particles in the loop - the 
last term of ()1.1|) . This generalises the known special cases with adjacent negative-helicity 
gluons, and the five-gluon non- adjacent amplitude. Moreover, this is the first example of 
the application of the MHV diagram approach to non-supersymmetric loop amplitudes, 
and provides further evidence that all cut-constructible (parts of) amplitudes may be 
derived using standard MHV diagrams. Of course, it would be extremely interesting to 
extend the MHV diagram method to obtain the rational pieces. This might require the 
construction of suitable MHV vertices where the off-shell legs are continued to 4 — 2e 
dimensions or the inclusion of additional effective vertices as proposed in [5]. 

^The paper [5] discusses the twistor structure of some non-supersymmetric one- loop amplitudes and 
the possible role of additional vertices in these models. A recent paper [18] has also developed a gener- 
alised MHV diagram construction for scattering amplitudes involving a Higgs boson and gluons. These 
amplitudes are described in terms of a tree-level, non-supersymmetric effective interaction which arises 
by integrating out a heavy top quark in one-loop diagrams. 
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The plan for the rest of the paper is as follows. In Section 2 we present the formal ex- 
pression for the one-loop MHV diagrams with a complex scalar running in the loop, which 
we use in Section 3 to rederive the known amplitude when the negative-helicity gluons are 
in adjacent positions. In Section 4 we derive a compact expression for the amplitude in 
the case where the negative-helicity gluons are in arbitrary positions. Our final result is 
given by Eq. ()4.22|1 . We also briefly comment on the twistor space structure of our result. 
Section 5 is devoted to some consistency checks of our general amplitude. Specifically, we 
show that it correctly incorporates the adjacent case result [19], also directly reproduced 
in Section 3, and the cut-constructible part of the five-gluon amplitude computed in [20]. 
Finally, we show that our expression has the expected infrared singularities. 



For further related work on the string theory side, and on the gauge theory side, 
see [22]-[27] and [28]-[34] respectively. 



2 The scalar amplitude 



In complete similarity with the A/' = 4 and A/'=l cases, see e.g. [15], we can immediately 
write down the expression for the scalar amplitude in terms of MHV vertices as 



^r'^^ = Yl ldMAi-lf,m,,...,i-,...,m2,lt] 

mi, 7712, ± 

■Ai-q,m2 + l,...,r,...,m, - 1,/f) , (2.1) 

where the ranges of summation of mi and 1712 are 

j + l<mi<i , i<m2<j - I . (2.2) 

The typical MHV diagram contributing to A^^^^^"^, for fixed mi and 1712, is depicted in 
Figure 1. The off-shell vertices A in p.ljl correspond to having complex scalars running 
in the loop. It follows that there are two possible helicity assignments^ for the scalar 
particles in the loop which have to be summed over. These two possibilities are denoted 
by ± in ()2.1|) and in the internal lines in Figure 1. It turns out that each of them gives 
rise to the same integrand for ()2.H) . 

_ . .tree {^2 ^2 + 1) K-l 7^) {t hf {j hf {l {j h? , . 

" {^J)' {mih) {mi-lh) {m2l2) {m2 + ll2) {hky ' ^ ' ' 

A crucial ingredient in ()2.H1 is the integration measure dAi. This measure was constructed 
in [6], using the decomposition L := I + zt] for a non-null four- vector L in terms of a null 
vector / and a real parameter z. tj is a null reference vector, which disappears in the final 

^For scalar fields, the "helicity" simply distinguishes particles from antiparticles (see, for example, [21]). 
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Figure 1: A one-loop MHV diagram with a complex scalar running in the loop, computed 
in Eq. ^2.\\ . We have indicated the possible helicity assignments for the scalar particle. 



result. We refer the reader to Sections 3 and 4 of [6] for the construction of this measure 
(also reviewed in Section 3 of [15]), and here we merely quote the result: 

dz 

dM = — rf'-'^LIPS(/2,-/i;PL;.) , (2.4) 
z 

where := U + Zifj, i = 1,2 and z := zi — Z2. Thus the integration measure dj^ 
decomposes into the product of a Lorentz-invariant two-particle phase space measure and 
a dispersive measure dz/z. The momentum Pl-z flowing in the phase space measure is 

Pl;z := Pl- z^l . (2.5) 

The interpretation of dz/z as a dispersive measure follows at once when one observes 
that [6] 



z Ph. - PI 



(2.6) 



In order to calculate ()2.H) . we will first integrate the expression ()2.3|) over the Lorentz 
invariant phase space (appropriately regularised to 4 — 2e dimensions), and then perform 
the dispersion integral. For the sake of clarity, we will separate the analysis into two 
parts. Firstly, we will present the (simpler) calculation of the amplitude in the case where 
the two negative-helicity gluons are adjacent. This particular amplitude has already been 
computed by Bern, Dixon, Dunbar and Kosower in [19] using the cut-constructibility 
approach; the result we will derive here will be in precise agreement with the result in 
that approach. Then, in Section 4 we will move on to address the general case, deriving 
new results. 
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3 The scattering amplitude with adjacent negative- 
helicity gluons 



The adjacent case corresponds to choosing i = mi, j = mi — 1 in Figure 1. Therefore we 
now have a single sum over MHV diagrams, corresponding to the possible choices of 
We will also set i = 2, j = 1 for the sake of definiteness, and m2 = m. 



After conversion into traces, the integrand of ()2.H) takes on the form: 



2Hh ■ k2)^h ■ [ ih-m+l) ik-m) 

where we note that [li ■ I2) = —Pj^.^/2 by momentum conservation. 



(3.1) 



The next step consists of performing the Passarino-Veltman reduction [35] of the 
Lorentz invariant phase space integral of 1)3.11) . This requires the calculation of the three- 
index tensor integral 

X^^^(m,Pi;,) = LlIPS(/2,-/i;Pl;.)|^ . (3.2) 

This calculation is performed in Appendix A. The result of this procedure gives the 
following term at 0{e^), which we will later integrate with the dispersive measure: 

fscalar ^2 i^^+ih h ^rn f^z)? i tl+ih h f^z , 2(A;i ■ A;2) 



^scalar _ _ (p^^ \-<^ t^^-hvr-i- rm 7 I "^-t-vri 7 i^,z rmj _^ 



— (m m + 1) , (3.3) 

and we have dropped a factor of AnXA^'^^^ on the right hand side of ()3.3)1 . where A is 
defined in (jdl)) . We can reinstate this factor at the end of the calculation. We also 
notice that ()3.3j) is a finite expression, i.e. it is free of infrared poles. 

An important remark is in order here. On general grounds, the result of a phase space 
integral in, say, the P^-channel, is of the form 

Ae) = {-P'r ■ m , (3.4) 

where 



/(e) = ^ + /o + /le + ■■■ , (3.5) 

and fi are rational coefficients. In the case at hand, infrared poles generated by the phase 
space integrals cancel completely, so that we can in practice replace ()3.5|) by /(e) 
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/o + /ic + ■ ■ ■ • The amplitude A is then obtained by performing a dispersion integral, 
which converts ()3.4|) into an expression of the form 



A{e) 





9o 



g,\og{-P') + g, + 0{e) , 



(3.6) 



e 



e 



where g{e) = go + fi'ic + ■ ■ ■ , and the coefficients gi are rational functions, i.e. they are free 
of cuts. Importantly, errors can be generated in the evaluation of phase space integrals 
if one contracts (4 — 2e)-dimensional vectors with ordinary four-vectors. This does not 
affect the evaluation of the coefficient go := g{e = 0), and hence the part of the amplitude 
containing cuts is reliably computed; but the the coefficients gi for i > 1, in particular 
gi, are in general affected. This implies that rational contributions to the scattering 
amplitude cannot be detected [19] in this construction. A notable exception to this is 
provided by the phase space integrals which appear in supersymmetric theories. These 
are "four-dimensional cut-constructible" [19], in the sense that the rational parts are 
unambiguously linked to the discontinuities across cuts, and can therefore be uniquely 
determined.^ This occurs, for example, in the calculation of the A/' = 4 MHV amplitudes at 
one loop performed in [6] . In the present case, however, the relevant phase space integrals 
violate the cut-constructibility criteria given in [19]^, since we encounter tensor triangles 
with up to three loop momenta in the numerator. Hence, we will be able to compute the 
part of the amplitude containing cuts, but not the rational terms. In practice, this means 
that we will compute all phase space integrals up to C(e°), and discard 0(e) contributions, 
which would generate rational terms that cannot be determined correctly. 

After this digression, we now move on to the dispersion integration. In the center of 
mass frame, where P^-z '■= Pl;z{^,0), all the dependence on P^-z in fl3.3|) cancels out, as 
there are equal powers of Pl-z in the numerator as in the denominator of any term. As a 
consequence, the dependence on the arbitrary reference vector r] disappears (see [16] for 
the application of this argument to the A/" = 1 case). Using ()2.6p in order to re-express 
dz/z in terms of the relevant dispersive measure, we see that we are left with dispersion 
integrals of the form 



The momentum flow can be conveniently represented as in Figure 2, where we define 



r fin' 1 
I{Pl) := ^^,{sT = -[vrecsc(7re)](-P|)-^. (3.7) 






P '■— <]2,m~l 7 Q '■— Q'm+1,1 — — (l2,m j 



(3.9) 



^For more details about cut-constructibility, see the detailed analysis in Sections 3-5 of [19]. 
"^An example of an integral violating the power-counting criterion of [19] is provided by I) A. 3(1 . 
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and 



P2 



Ef=pi We also have Pl : = 



92. 



-g. 



Now we wish to combine the terms in the first hne of ()3.8|) with those in the second 
line. Since flH.8|l is summed over m, we simply shift m + 1 ^ m in the terms of the second 
line. Let us now focus our attention on the second term in (similar manipulations 
will be applied to the first term). Writing the m ^ m + 1 term explicitly, we obtain a 
contribution proportional to 



[tr+(^i M^h) 



[tr+(^i^2 f-m+ih)\ 



By shifting m+ 1 ^ m in the second term of (jH.lOj) . we convert its Pl to Pl q2,m-i 
(whereas, in the non-shifted term, Pl = —Q)- The expression ()3.10|) then reads 



(3.10) 



P 



-Q 



2\-t 



where we used tr+(, 

m-g = -(i/2)(g2 



^1 11^2 



(m-g)2 

)) = -ti+{^i ^2 ^mf) and g 



P 



2\-e 



m 



(3.11) 

-P ■ m. Notice also that 



m 




Figure 2: A triangle function contributing to the amplitude in the case of adjacent negative 
helicity gluons. Here we have defined P := qj,m-i, Q '■= <im+i,i = ~1j,m (in the text we set 
i = 1, j = 2 for definiteness). 

Next we re-instate the antisymmetry of the amplitudes under the exchange of the 
indices 1^2 (which is manifest from equation ()2.3|) ). Doing this we get 

1 



(3.12) 



2{ki-k2){m-Q) tr+(^i^2^. 



tr+(^i ^2 
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Following similar steps for the first term in ()3.8j) . we arrive at the following expression for 
the amplitude before taking the e ^ limit: 



2,e , 



(3.13) 



where 



^tree 


TT 




6 ' 


^tree 


TT 




' 3 ' 



tr+(|^i 1^2 |^m^m,l) - tr+(|^i |^2^m,l |l^m) (m, g2,m-l , g2,m) , 
[tr+(|^l ^2 im,l)\ ^tr+(|^i 1^2 ^m,l ^ni) - 

) , (3.14) 

and ti follows from the definitions below equation ()4.8j) . In order to write ()3.14j) in a 
compact from, we have introduced e-dependent triangle functions [15] 



e (Q2 - P2)r 



(3.15) 



where p + P + Q = 0, and r is a positive integer.^ 

We can now take the e — limit. As long as and are non- vanishing, one has 
liniTW(p, P, Q) = T^^\p, P,Q), P V , Q V , (3.16) 

where the e-independent triangle functions are defined by 

-)2 / r>2\ 



IP,^,<VJ • (Q2_p2)r- 



(3.17) 



If either of the invariants vanishes, the limit of the e-dependent triangle gives rise to 
an infrared- divergent term (which we call a "degenerate" triangle - this is one with two 
massless legs). For example, if = 0, one has 



T.(p,P,Q)|q.=o 



1 (_p2)-e 
7 P2 







(3.18) 



The two possible configurations which give rise to infrared divergent contributions corre- 
spond to the following two possibilities: 



a. ^2,™-! = ^2 (hence q2m-i = 0)- this case we also have qim = ^ 



[2] 



-?2. 



ki (hence g'im = 0)- Therefore = tn 



[2] 



^For r — 1 we will omit the superscript (1) in T^^^ 
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We notice that infrared poles will appear only in terms corresponding to the triangle 
function T. Indeed, whenever one of the kinematical invariants contained in T^^^ vanishes, 
the combination of traces multiplying this function in ()3.14p vanishes as well. 

In conclusion, we arrive at the following result, where we have explicitly separated out 
the infrared-divergent terms:^ 

AT'^' = Aoics + A, + A2, (3.19) 

where 

1 1 r r„, 1 

(3.20) 



1 1 

Ai = ^[tr+(^i^2^m^m,l) - tr+(|^l|^2^m,l^m) T{m, q2^m~l, q2,m) 



tl m=4 



n-1 



A2 = i^*■■'=<=-^^[[tr+(^l^2^„^n^,l)]V(l^l|^2^„^,l^r 



m=A 



- tr+(|^i|^2|^m^m,l)[tr+(|^l|(^2^m,l|^m)]^ ^^^^ (m, g2,m-l , ^2,™) • 

More compactly, we can recognise that ^poics and Ai reconstruct the contribution of an 
J\f=l chiral supermultiplet, and rewrite ()3.19|1 as 

1 11""^ 

^scalar ^ _^^=l,ch.al ^ -^f- ^ iJ™ T(3) (m, g2,^_l , ^2,..) , (3.21) 



where 



m=4 



BT2 = [iT+{hMn.im,i)\ tr+(^i^2^™,i W (3.22) 

- [tr+(|^i ^2 im,l ^m)] ^tr+(|^i 1^2 im,l) ■ 



and 



1 1 " 

^^=l,chiral ^ _jtvee_ ^ |^tr+ (|^i |^2 ^m,l) " tr+ (|^i |^2 ^m,l |^m) T{m,q2,m-l,q2,m) ■ 

^1 m=3 

(3.23) 

This is our result for the cut-constructible part of the n-gluon MHV scattering amplitude 
with adjacent negative-helicity gluons in positions 1 and 2. This expression was first 
derived by Bern, Dixon, Dunbar and Kosower in [19], and our result agrees precisely with 

factor of —AttX will be understood on the right hand sides of Eqs. (|3.19(l . H3.21|l . I|3.23|l . where A is 
defined in l|C.l|l . 
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this. A remark is in order here. In [19], the final resuh is expressed in terms of a function 



L2{x) 



log X — {x — l/x)/2 



(3.24) 




which contains a rational part — (x — l/a:)/2(l — x)^. This rational part removes a spu- 
rious third order pole from the amplitude, but with our approach we did not expect to 
detect rational terms in the scattering amplitude, and indeed we do not find such terms^. 
Furthermore, we do not find the other rational terms which are known to be present in 
the one-loop scattering amplitude [20]. 



4 The scattering amplitude in the general case 

The situation where the negative- helicity gluons are not adjacent is technically more 
challenging. Our starting point will be ()2.3|1 . to which we will apply the Schouten identity 
(see Appendix D for a collection of spinor identities used in this paper). Eq. ()2.3|) can 
then be written as a sum of four terms:^ 

C(mi,m2 + 1) - C{mi,m2) - C(mi-l,m2 + l) + C{mi - 1,1712) , (4.1) 



The calculation of the phase space integral of this expression is discussed in Appendix B. 
The result is 



where 



C{a,b) 



(4.2) 




1 ti+{^iH) 
3 (a ■ 6) 



tr+(/ifi;,^) 



2\tT^MHh^ 



+ 



(a ^ b) 




1 il+{iiH)il+{^iif>) i^L.^if 

2 (a ■ bf {Pl;. ■ af 



+ {a^b) 



tr+(/^i#)tr+(/^i/^^) tr+(/^i^/^)tr+(/if^.,^) 
(a ■ bf {Pl;, ■ a) 



+ (a 




''in our notation L2 corresponds to T^^\ which, however, lacks a rational term. 

®We drop the factor of — i.A^''°° from now on and reinstate it at the end of the calculation. 
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+ 



"tr_i 



(a • by 



log 1 - 



b) 



N 



p2 



(4.3) 



where := N{P) := (a • b)P^ — 2(P ■ a){P ■ b), and we have suppressed a factor of 
-47rA(-P|)~' ■ [2^{i ■ j'f]-^ on the right hand side of (jOI), where A is defined in ^T\\ . 
We notice that ()4.3|) is symmetric under the simultaneous exchange of i with j and a with 
b. This symmetry is manifest in the coefficient multiplying the logarithm - the last term 
in ()4.3|) : for the remaining terms, nontrivial gamma matrix identities are required. For 
instance, consider the terms in the second line of ()4.3|1 . These terms are present in the 
adjacent gluon case ()3.3j) . and it is therefore natural to expect that the trace structure of 
this term is separately invariant when ^ j and a ^ b. Indeed this is the case, thanks 
to the identity 



32(z-j) 



tr. 



+ tr. 



tr. 



L;z) 



{Pl;z 



L-z) 



tr. 



■ a) 



(4.4) 



{Pl;z ■ 



+ 



Similar identities show that the third and fourth line of (H. 
simultaneous exchange i ^ j and a ^ b. 



2(^-j) 
{Pl;z ■ a) 

2(^-j) 
{Pl;z ■ a) 

are invariant under the 



The next step is to perform the dispersion integral of ()4.3p . i.e. the integral over the 
variable z. This appears in the terms involving P^.^ in ()4.3p . and in an overall factor 
{Pl z)~^ arising from the dimensionally regulated measure. 

The integral over the term involving the logarithm has been evaluated in [6] , with the 
result 

Notice that these terms were not present in the adjacent negative-gluon case considered 
in Section 3. 



Next we move on to the remaining terms in ()4.3|l . Inspecting their ^-dependence, we 
see that, in complete similarity with the adjacent case of Section 3, in each term there 
are the same powers of Pl-z in the numerator as in the denominator. Hence, in the centre 
of mass frame in which Pl-^ '■= Pl;z{^i^)i one finds that Pl-z cancels completely. Note 
that this also immediately resolves the question of gauge invariance for these terms - this 
occurs only through the rj dependence in Pl-^z = Pl — zrj. Furthermore, the box functions 
coming from ()4.5p are separately gauge invariant [6] . The conclusion is that our expression 
for the amplitude below, built from sums over MHV diagrams of the dispersion integral 
of ()4.3|) . will be gauge invariant. Moreover, apart from ()4.5|) . the only other dispersion 
integral we will need is that computed in ()3.7|1 . 
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It follows from this discussion that the result of the dispersion integral of ()4.Hj) is 
(suppressing a factor of — 47rA(— Pl)^" ■ [2^(i ■ j)^]^^ ■ [vre csc(7re)]): 

— /rf4-2^LIPS(/2,-/i;PL;.) C{a,b) 



-pir 



1 tr_ 
3 ~ 



(a-b) 



tr. 



(Pl ■ a)3 {Pl ■ a) 



— {a ^ b) 



1 tr+(/i/^^)tr+( 

2 (a ■ by 



(Pl ■ ay 



+ {a^b) 



tr. 



tr. 



(a ■ by 



{Pl ■ a) 



b) 



tr. 



"trj 



(a ■ by 



Lis 1 - 



N{Pl 



-Pr 



(4.6) 



Now, due to the four terms in ()4.H1 . the sum over MHV diagrams will include a signed 
sum over four expressions like ()4.6|) . Let us begin by considering the last line of ()4.6|) . 
This is a term familiar from [6] and [15], corresponding to one of the four dilogarithms in 
the novel expression found in [6] for the finite part 5 of a scalar box function. 



B{s,t,P\Q') = Li2 ( 1 - ]^P^J + Lis ^1 - 



N{P) J \ N{P) 

with s := (P + ay, t := {P + by, and P + Q + a + b = 0. By taking into account the four 
terms in ()4.H) and summing over MHV diagrams as specified in ()2.H) and ()2.2j) . one sees 
that each of the four terms in any finite box function B appears exactly once, in complete 
similarity with [6] and [15], so that the final contribution of this term will be^ 



1 2 

mi,m2 — 1' 9'mi+l,m2' 9'mi+l,m2 — 1' ^m2+l,ra\ — i) , (4.8) 

mi=j+l m2=i+l 

where ^ := (pj + pj+i + ■ ■ ■ + Pi+fc_i)^ for A; > 0, and = tl""*^^ for /c < 0. In writing 
()4.8j) . we have taken into account that the dilogarithm in ()4.(jj) is multiplied by a coefficient 
proportional to the square of b^l^^^, where 

[{ki + kjyy [{kmi + kra2y 

^We multiply our final results by a factor of 2, which takes into account the two possible helicity 
assignments for the scalars in the loop. 



"mima ^ [rt. _L .^212 r^i. i ^212 ' ^^'^ ' 
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We notice that femim2 



is the coefficient of the box functions in the one- loop A/" = 1 MHV 
amphtude, originally calculated by Bern, Dixon, Dunbar and Kosower in [19], and derived 
in [15, 16] using the MHV diagram approach for loops proposed in [6]. Furthermore, 



rriH-l 



m, 



,+1 ^ 



mr,-1 



Figure 3: A box function contributing to the amplitude in the general case. The negative- 
helicity gluons, i and j , cannot be in adjacent positions, as the figure shows. 

we observe that ^^^^2 holomorphic in the spinor variables, and as such has simple 
localisation properties in twistor space. Indeed, from ()4.9p it follows that 

,ij ^ {inii) (ims) {j mi) (jms) 

'"^'"^ (ij)2(mim2)2 ■ ^ ' 

Summing over the four terms for the remainder of ()4.6|) can be done in complete similarity 
with Section 4 of [15].^° We will skip the details of this derivation, and will now present 
our result. 



In order to do this, we find it convenient to define the following expressions: 
{ijm2 + lmi) (ij 777,2 mi) 



m\m2 



((7772 + 1) • rrii) 
-2 [ij] {nil i){mij) 



(7772 ■ rrii) 

(7772 7772 + 1) 

{m2 + 1 777i) {nil ^2) 



(4.11) 



^"^In Section 3 we have illustrated in detail how this sum is performed for the simpler case of adjacent 
negative-helicity gluons. 
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m\m2 



TV 

in\m2 



[ij mi 1712 + ^2 + 1 ^i) {ij nT-i ^2){ij ^2 ^i) 



((m2 + 1) ■ miY 

{ij mi m2 + "^2 + 1 ^i) 

((m2 + 1) • miY 



(m2 ■ mi)2 ' 

{ij mi m2Y{ij m2 mi) 
{m2 ■ miY 



(4.12) 



, (4.13) 



where for notational simplicity we set {ai 02 0304) := tT^{^i^2i34i) the above. We also 
note the symmetry properties 



m\m2 



m\m2 ' 



Q3^ 

m\m2 



m\m2 



(4.14) 



The momentum flow is best described using the triangle diagram in Figure 4, where we 
use the following definitions: 



P ■ C[m2+l,mi — l 
Q • Qmi+l,m2 • 



m2 1 



(4.15) 



The triangle in Figure 5 also appears in the calculation, and can be converted into a 
triangle as in Figure 4 - but with i and j swapped - if one shifts mi — 1 — mi, and then 

swaps mi <-> 772,2. 



We then introduce the coefficients 



'^mim2 


■= 2-^ 




j) 


-4 Aij 

m\m2 


{ijmiQf{i]Qmi) - (- 


:jmi Q){ijQmif 


,(4.16) 


'^m\m2 


:= 2-^ 




j) 


-4 Aij 

m\m2 


{ijmiQf - {ijQmif 


5 


(4.17) 




= 2-« 


{^ 


j) 


-4 qij 


{ijmiQf + {ijQmif 


5 


(4.18) 


jij 

m\m2 


:= 2-« 


{^ 


j) 


-4 

mim2 


{ijQmi) + P^^^^{ijmi 


Q) ■ 


(4.19) 



We will also make use of the e-dependent triangle functions introduced in ()3.15|) . whose 
e — > limits have been considered in ()3.16|) - ()3.18|) . This is in order to write a compact 
expression which incorporates also the infrared- divergent terms. 

^^The infrared-divergent terms will be described below, and used to check that our result has the 
correct infrared pole structure. 
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mim2 +1,7712' ^mi+l,m2 — 1' ^m2+l,m\ — li 



We can now present our result for the one- loop MHV amplitude ()2.H1 ^^: 

j-i i-i 

2 

mi=j+l m2=i+l 
mi=j+l m2=i 
mi=j+l m2=i 



(4.20) 



where on the right hand side of ()4.2()j) a factor of —AttX is understood, where A is defined 
in ()C.1|) . We can also introduce the coefficient 



{ijm2 + lmi) {ij 7712 mi) 



_((m2 + l)-mi) (m2-mi)J [(« + jT] 



{ijmiQ) - jijQmi) 

\2]2 



(4.21) 



which already appears as the coefficient multiplying the triangle function T in the A/'=l 
amplitude, (see e.g. Eq. (2.19) of [15]), and rewrite ()4.2n|l as 



mi,m2 — 1' '?mi+l,m2' +I,m2 — 1' ?m2+l,mi — l) 

^ mi=j+l m2=j+l 



scalar 



mi=j+l m2=i 



(^.■rn,g)(.,Qr.,) ^,3)(^^^p^^) ^ T(m„P,Q) 



2(^-jT 



+ 2 5^ 5^ [5^^^^T(2)(mi,P,Q) + J^^^^T(mi,P,Q) 

mi=i+l m2=i 



+ {i ' ^ j] 



(4.22) 



Several remarks are in order. 



1. As usual, the variables ^2-1' 1rni+im2 correspond to the s- and t-channel of the 
finite part of the "easy two-mass" box function with massless legs mi and m2, and 
massive legs ql^+i ^^^_i (Figure 3). 

^^We thank Lance Dixon for pointing out some typos, both in these equations and elsewhere, in an 
earher version of the paper. 



15 



^1 




Figure 4: One type of triangle function contributing to the amplitude in the general case, 
where i & Q, and j G P. 

2. Compared to the range for mi and m2 indicated in ()2.2|1 . we have omitted mi = i 
in the summation of the triangles, as for this value the coefficients A, S, I defined 
in ()4.16p - ()4.19|) vanish. Notice also that we have i E Q and j G P. 

3. In the case of adjacent negative- helicity gluons, the only surviving terms are those 
containing the coefficient d^^^^, on the second line of ()4.2m) or ()4.22|) . We will 
return to this point in Section 5. 

4. We comment that, in contrast to the adjacent case (see ()3.21|) ). in the general case 
the N = 1 chiral amplitude does not separate out naturally in the final result, as 
one can see from the coefficient of the box function B in ()4.20|) . 

Next we wish to separate explicitly the infrared divergences from ()4.2Up . We can 
immediately anticipate that there will be four infrared-divergent terms, corresponding 
to the four possible degenerate triangles. Two of these degenerate triangles occur when 
either P^ or happen to vanish. The other two originate from the i ^ j swapped terms. 

Let us consider first the terms arising from the summation with i ^ j unswapped. 
When = 0, it follows that mi = i — 1 and m2 = i (see Figure 4). When P^ = 0, it 
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Figure 5: Another type of triangle function contributing to the amplitude in the general 
case. By first shifting mi — 1 — » mi, and then swapping mi ^ m2, we convert this into 
a triangle function as in Figure 5 - hut with i and j swapped. These are the triangle 
functions responsible for the i ^ j swapped terms in ()4.2()|1 - or ()4.22|1 . 



follows that mi = j + 1 and m2 = j — 1 (see Figure 5). Hence 



TM(p,p,g) 

T«(p,P,Q) 



1 (-^. 



(41)^ 



(4.23) 



The infrared-divergent terms coming from Q = are then easily extracted, and are 



1 

2^ 



(-i?ii)-^4(z.j) 



(iji — li + l) 



{i i i + \ i — 1) 



(4.24) 



((^ + 1) ■ (^ - 1)) 



{iji + li - l){iji- l'i + 1) 
((^ + l)-(^-l))' 



and from = 



^ (-tf^)-M(.-j) 



2e 



{ijj - 1 j + 1) 
((j + 1) ■ (j - 1)) 
{ijj + Ij - 1) 



(4.25) 



((j + 1) • ij - 1)) 



(^ ■ j) + 



(O' j + Ij - 1)(0' j - 1 j + 1) 
((J + 1)-(J-1))' 
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Likewise, from the "swapped" degenerate triangles we obtain the following infrared- 
divergent terms: 



1 

2^ 



{ijj + 1 j - 1) 
((j + 1) ■ (j - 1)) 



(4.26) 



8. s9 (ijj-lj + 1) , , (z 7 7 - 1 7 + ?■ 7 + 1 7 - 11 



and 



1 

Ye 



-t 



{iji + li — l) 



(4.27) 



8 .^2 o (ijz-li+1) ^. (zjz-H + l)(zji + li-l) 



4.1 Comments on twistor space interpretation 



We would like to make some brief comments on the interpretation in twistor space of our 

resuh srm . 



1. As noticed earlier, the coefficient b^l^^^ appears already in the Af=l chiral super- 
multiplet contribution to a one-loop MHV amplitude, where it multiplies the box 
function. It was noticed in Section 4 of [5] that b^l^^^ is a holomorphic function, 
hence it does not affect the twistor space localisation of the finite box function. 

2. The coefficient c^^^^ ^^^^ appears in the Af=l amplitude, as the coefficient of the 
triangles (see e.g. Eq. (2.19) of [15]). Its twistor space interpretation was considered 
in Section 4 of [5], where it was found that c^^m2 support on two lines in twistor 
space. Furthermore, it was also found that the corresponding term in the amplitude 
has a derivative of a delta function support on coplanar configurations. 

3. The combination c^^^^ Q) (O Q"^i)/(^ " j)^ already appears in the case of 
adjacent negative- helicity gluons. The localisation properties of the corresponding 
term in the amplitude were considered in Section 5.3 of [5], and found to have, 
similarly to the previous case, derivative of a delta function support on coplanar 
configurations. 

4. On general grounds, we can argue that the remaining terms in the amplitude have a 
twistor space interpretation which is similar to that of the terms already considered. 
The gluons whose momenta sum to P are contained on a line; likewise, the gluons 
whose momenta sum to Q localise on another line. 

^We thank Dave Dunbar for discussions on this point. 
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We observe that the rational parts of the amphtude are not generated from the MHV 
diagram construction presented here. Such rational terms were not present for the M = 1 
and A/" = 4 amplitudes derived in [6,15,16]. However, for the amplitude studied here, 
rational terms are required to ensure correct factorisation properties [19]. 



5 Checks of the general result 

In this Section we present three consistency checks that we have performed for the result 
()4.20j) (or ()4.22|) ) for the one-loop scalar contribution to the MHV scattering amplitude. 
These checks are: 

1. For adjacent negative- helicity gluons, the general expression ()4.2()|1 should reproduce 
the previously calculated form ()3.2H) . 

2. In the case of five gluons in the configuration (l~2+3~4+5^), the result ()4.20|) should 
reproduce the known amplitude given in [20]. 

3. The result ()4.20p should have the correct infrared-pole structure. 
We next discuss these requirements in turn. 



5.1 Adjacent case 

The amplitude where the two negative- helicity external gluons are adjacent is given in 
Section 7 of [19] and was explicitly rederived in Section 3 of this paper by combining 
MHV vertices, see Eq. ()3.2H) . It is easy to show that our general result ()4.22|) reproduce 
correctly ()3.21|1 as a special case. 

To start with, recall that our result ()4.22j) is expressed in terms of box-functions and 
triangle functions, see Figure 3 and Figures 4, 5 respectively. In the adjacent case, the 
box functions are not present. Indeed, in the sum ()4.8|) the negative- helicity gluons can 
never be in adjacent positions (see Figure 3). 

Next, we focus on the triangles of Figure 4. In terms of these triangles, requiring i and 
j to be adjacent eliminates the sum over m2, as we must have m2 = i and m2 + 1 = j. 
Moreover, in this case Q = qmi+i,i, P = Q'j.mi-i and one has 
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Qmim2 







Tmim2 







(5.1) 



for 777.2 = i, and m2 + 1 = j. Similar simplifications occur for the swapped triangle. Hence 
the only surviving terms are those in the second line of ()4.20j) (or ()4.22j) ). and it is then 
easy to see that they generate the same amplitude ()3.8|) already calculated in Section 3. 



5.2 Five-gluon amplitude 

The other special case is the non-adjacent five-gluon amplitude (l~2+3^4+5^), given in 
Eq. (9) of [20]. This amplitude may be written as cr^Atree times^^ 



^-gl0g(-S34) + 2^(2.5)4(1.3)4 ^(^51,^12,0,^34) 



1 tr+(/^;^^) tr+(/^^^) 
3 24(2 ■5)(1- 3)4 



tr4/ ^ f> MWfsjl ^ ^^^(^ ^ ^ k)g(W£5il 



(Si2 — S34) 



+ 



1 



1 



3 23(1-3)2 



tr+(/^^^)tr+(/^^. 



^2 log(g34/g5l) 
(S34 - S5i)3 



+ 



tr+(/^^^)nr+(/^^^) 
26(2 -5)2(1 -3)4 

tr+(;^;^^)2 tr+(/^^^) 
26(2-5)3(1-3)4 



I0g(si2/S34) I0g(s34/S5i) 



(Sl2 - 534)^ (S34 - S5l)^ 
I0g(si2/S34) , I0g(s34/S5i) 



3 22(1-3) 



(Sl2 - S34, 

log(g34/g5l) 
(S34 - S51) 



IS34 - S51, 



+ (1,4)^(3,5) 



(5.2) 



where the interchange on the last line applies to all terms above it in this equation, 
including the first two terms, and the box function B is defined in ()4.7|1 . In deriving this 
from [20], we have used the dilogarithm identity 

fl — r\ fl — s\ fl — sl — r 
Li2(l - r) + Li2(l -s) + log(r) log(s) = Li2 ( j + Li2 ( j - Li2 ( — — 

(5.3) 



-'^^The derivation in [20] used string-based methods, which affects the coefficient of the pole term. In 
(|5.2|l we have written the pole coefficient which matches the adjacent case. 
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We have checked exphcitly that our expression for the ra-gluon non- adjacent amphtude 
()4.20|) . when speciahsed to the case with five gluons in the configuration (1^2+3^4+5^), 
yields precisely the resuh fl5.2|) above. For the terms involving dilogarithms, this is easily 
done. For the remaining terms, which contain logarithms, a more involved calculation 
is necessary using various spinor identities from Appendix D. A straightforward method 
of doing this calculation begins with the explicit sum over MHV diagrams in this case, 
isolating the coefficients of each logarithmic function such as e.g. log(si2), and then check- 
ing that these coefficients match those in ()5.2|) . The remaining 1/e term arises from the 
following discussion. 



5.3 Infrared-pole structure 



The infrared-divergent terms (poles in 1/e) can easily be extracted from ()4.24j) - ()4.27j) by 

[91 ______ 

simply replacing {—tr )^^ — > 1 (r = z — — l,j). Consider first the terms in (I4.25|l 

and ()4.2f)j) . After a little algebra, and using 

(zjj + lj-1) + (zjj-H-1) = 4(z-j)((j-l)-(j + l)) , (5.4) 

one finds that these two contributions add up to 

64 

-^i^-jf- (5.5) 
3 e 

Similarly, the pole contribution arising from ()4.24j) and ()4.27j) gives an additional contri- 
bution of — (64/3 e) {i ■ j)^. Reinstating a factor of —2 • 2^^{i ■ ■ Atree, we see that the 
pole part of ()4.20|) is simply given by 

-^scalar I £_polc ~ ^ • (5-6) 

Hence our result ()4.2()|1 has the expected infrared- singular behaviour. 
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Appendix A: Passarino-Veltman reduction 



In Section 2 we saw that a typical term in the cut-constructible part of the Yang-Mills 
amplitude is the dispersion integral of the following phase space integral: 



tr+(^l ^2 h; zf 2) tl+ih hhfL;z) tl+ {h h h) 



The goal of this Appendix is to perform the Passarino-Veltman reduction [35] of (jA.lj) . 
To this end, we rewrite C(m) as 

- (A;i . hriPlJ' ^ irn,F,;z), (A.2) 



where^^ 



/M 11/ JP 

'2 '2 '■2 



J [I2 ■ m) 

On general grounds, X'^^P{m, Pl) can be decomposed as 

+ {m^PlPi + Pi:m'^P^ + Pi:P^W)J3 + P^PIPIT, 

+ {r]>"'mP + r^^^m'^ + r/^'^m^) J5 + (r^^'^Pf + r/^^Pj^ + r/^^P^) Jg , (A.4) 

for some coefficients Tj,i = 0, ...,6. One can then contract with different combinations of 
the independent momenta in order to solve for the Xj. Introducing the quantities 



A 
B 
C 
D 
E 
F 



m^m^PLpI'"'' , 
m.PL.PLpI^'"' , 

il,,mpl^^p = 0, 

VpuPLpI^"' = 0, (A.5) 



the result for the Passarino-Veltman reduction of {Ii, . . . in the basis {A, . . . , D} is: 

J2 = {5{Plf/{2{m.PLf),-6Pl/{m.PL)\3/{m.PLf,0), 

J3 = (-2P|/(m-Pi)\3/(m-Pi)3,0,0) , 

J4 = l/(m■Pi)^ 0,0,0) 

= -(P^)V(2(m■Pz.r),3P|/(2(m■Pi)3),-l/(m■Pz.)^0) , 

Je = (P|/(2(m-Pif),-l/(m-P,.)2,0,0) . (A.6) 
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For the rest of this Appendix we drop the subscript z in P^-z for the sake of brevity. 
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We omit the decomposition for Xi, as the corresponding term in (jA.4|l drops out of all 
future expressions due to fc^ = 0. 

Finally, using the methods of [15] and the results of Appendix C, the integrals in ()A.5|1 
are found to be, keeping only terms to order C(e°), 

A = (m-Pi)2^7rA, (A.7) 
B = Pl{m-PL)7rX, (A.8) 



C = {PlfnX, (A.9) 
(m ■ Pl) e 



D = -^7^^ -A, (A.IO) 



where 



7r2 



A := ^ . (A.ll) 

4i-r(i-e) 



Appendix B: Evaluating the integral of C(a, 6) 



The basic expression which arises in the MHV diagram construction in this paper is 

{^h){Jhy{^l2y{Jl2) {^a){Jb) 

We wish to integrate this expression over the Lorentz invariant phase space. We begin by 
simplifying it, using multiple applications of the Schouten identity. First note that using 
this identity twice, one deduces that 

{ab) = {ia){b]) + {ia){a])—— + {bj){ib)—— (B.2) 



{hl2){ab) 



+ {aj){tb) - {aj){tb) 



{ah){bl2) • 



Now use this identity in C{a,b). This generates five terms, which we will label (in 
correspondence with the ordering arising from the order of terms in (IB.2jl above) as 
Tj, z = 1, . . . , 4, and U. The Tj have dependence on the loop momenta such that we may 
use the phase space integrals of Appendix C to calculate them. The term U is more com- 
plicated; however, one may again use the identity ()B.2|) . generating another five terms, 
which we will label T5, . . . ,Ts, and V. Again, the expressions in Ti,i = 4, ...,8 may be 
calculated using the integrals of Appendix C. Finally, the term V may be simplified, here 
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using the identity ()B.2j) with i and j interchanged. This generates a further five terms, 
which we label Tg, . . . , T13. The explicit forms of these terms follow: 



and 



and 



rj. ^ tr+(^i^^)^tr+(^ifif2)tr+i/t/)>2n/ /-p o 

2^^{z-jy{a-by{k-l2)^iz-b)ia-h) ' ^ ' 

^ ^ tT+{^im^+i^iH)^^+i^ih h)^^^^ih ri)tr+(i ^^2^5) 



^ ^ tr+(^i/^^)^tr+vt/^/;;u + vt/^2n/ ^7 
28(z-j)^(a-6)3(/i./2) ' ^ ■ 

' 2^^{i-jY{a-hf{h-h){i-a){h-h) ' ^ • 



_ tr+(/i^/^)^tr^,.,^^, 



tro^fi i ^|5)2tr^fi i W)2tr+i^^ 1 



■ 11 



2^\i-jY{a-hY{i-a){h-h) 

ti^iiiihYtu 



(B.13 



- 2^{i.jY{a-hY ' ^^-^^ 

^ ^ tr+(^i;5^)^tr^(^i#)V(;5f2ri^) 

2io(^-jr(a-6)4(a-Zi)(6-/2) ' ^ ' 

The expression C(a, h) is then the sum of the terms Tj, i = 1, . . . , 13. 

Before performing the phase space integrals, it proves convenient to collect the result- 
ing expressions in pairs as Ti + T2, T3 + T4, T5 + Tg, T7 + Tg, Tg + Tu and Tio + T12, we 
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are led to the following decomposition: 

tr+(^ifif2)tr+(^if2ri)tr+(^iri^)tr+(ii/5f2) 



where 



Hi 



Ho 



7^4 



-C{a,b) 



2\i-j)\h-hnk-a){h-h) 



1 



(Hi + 



{h ■ h? {a ■ h) 

tr+(^i#)tr+(/i;6^)tr+(^ir2ri: 



n 



ih ■ k) (a ■ by 



ih-a) ik-b) 



ih ■ a) 



ik ■ b) 



(a ■ bf 



tr. 



tr. 



(h ■ a) 



ih ■ b) 



4(a ■ 6)4 (k ■ a) {h ■ b) 



(B.16) 



(B.17) 



Finally, we perform the phase space integrals of the above expressions, using the 
formulae in Appendix C below. One finds quickly that the divergent (as e — 0) part of 
the total expression is zero. The finite part, after further spinor manipulations, becomes 
the expression we have given in ()4.3|) . 

Appendix C: Phase space integrals 



The basic method which we use for evaluating Lorentz-invariant phase space integrals 
has been outlined in our earlier papers [6,15]. Here we will just quote the results which we 
need. In the following we will use a shorthand notation where / = /(i4~^''LIPS(/2, — /i; -Pl), 
and a common factor of 47rA(— Pl;^)^ is understood to multiply all expressions, where A 
is the ubiquitous factor 



TX2 

X := ^ . (C.l) 

4i-r(i-e) 

In the following we define a = {a ■ P), p = {b ■ P), N{P) = (a ■ b)P^ - 2(a ■ P)(6 ■ P) and 
drop the L; z subscripts on P^.^ for clarity. 



Firstly we quote the results from Appendix B of [15] up terms of order (9(e°): 

1 1 r I I 

ea ' 



(a ■ h) 



{b ■ h) e(3 



(C.2) 



/ 



1 



[a ■ li){b ■ k) 



N{P) 



i + L 
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where 

Prom this, we can derive recursively the following integrals (up to 0(e°)): 

= \p', 1 = - (C.3) 

j IK = J m = \ (p'P" - Iv'^P' 



f 11 

J (a-h 
J {b- 



p2 I p2 



) 2eQ;^ a a 



and 



3P^ 



/A* 7!^ p4 1 p2 

= -T^a'^a'^ + -^P'^P'^ + f^P^'^a'^) - ^ ,^ ^ 

(a ■ /i) 4ea3 2q; 2a^ Aa-^ Aa 



V''' , (C.4) 



It^lu p4 1 p2 'iP^ p2 

^2 = J—b''F-—P''P''-—PH''^ + ^b%'' + —n 



1X1/ 



^ = ±(2P>^-^a'^ + ^bA+'-^(p'^-^a'^ + ^b^ 

{a-h){b-l2) eN\ a (3 J N \ {a ■ b) {a ■ b) 

Finally, there are integrals involving cubic powers of loop momenta in the numerator. The 

first is 

/jt^wjP pi p2 1 p4 p2 

^AJ_ = p(Ma-aP) + pit^p-aP^ + —pt^p-pp _ —rj^^^a"^ - —rj^'^'pp^ /C.5) 
(a • h) 4q;3 4q;2 Sa^ Aa ^ ' 

where we have suppressed terms cubic in a as they prove not to contribute when this 
integral is contracted into the products of Dirac traces which appear in the expressions 
in Appendix B. The second cubic integral required is 

/JfJ-JvjP p4 p2 1 p4 p2 
^-AA = PH'^fe''^ H pit'p^b''^ -I pt^p^'pp r?^'*''^'') r^ip-^pp) fc 6~) 
{b-h) 4/33 ^ AI3^ ^3/3 8/32'' 4/3^ > 

again suppressing terms cubic in b which will not contribute. 

Appendix D: Spinor identities 



We collect here some formulae useful for the calculations presented in this paper. The 
Schouten identity is 

{ij){kl)^{ik){jl) + {il){kj) . (D.l) 
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Other identities are: 

[i3]{3i) = tr+(^,^,-) = 2(A;i-%) , (D.2) 

[i3]{3l)[lm]{mi) = ix+{}^4jUm) . (D.3) 

[ij]{jl)[lm]{mn)[np]{pi) = tY+{J^4jJf:4m}^n^p) , (D.4) 
for momenta /cj, kj, ki, km, kn, kp. We also have, for null momenta i, j, k, a, b, 

U-a) - (i-a) ■ 

For dealing with Dirac traces, we have the following identities'^ 

^^+{l^i^j^li^m) — 4(A;j • kj){ki ■ k^) — ^'^+{l^j^i^l^m) -I (^•'') 

tr+iijl^f)tr+i/ji^4) = 0, (D.8) 

tr+i^jlxf)tT+i^j4fi) = A{i-3)tr^i^j4f) . (D.9) 



""^^The appearance of a Greek letter such as inside a trace indicates that the relevant gamma matrix 
is to be inserted at that point. 
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